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Abstract 

Itai and Rodeh showed that, on the average, the communication of a leader election 
algorithm takes no more than LN bits, where L ~ 2.441716 and N denotes the size of the 
ring. We give a precise asymptotic analysis of the average number of rounds M(n) required 
by the algorithm, proving for example that M(oo) := lim M(n) = 2.441715879 . . ., where n 

n — >oo 

is the number of starting candidates in the election. Accurate asymptotic expressions of the 
second moment M' 2 ' (n) of the discrete random variable at hand, its probability distribution, 
and the generalization to all moments are given. Corresponding asymptotic expansions 
(n — > oo) are provided for sufficiently large j, where j counts the number of rounds. Our 
numerical results show that all computations perfectly fit the observed values. Finally, we 
investigate the generalization to probability t/n, where t is a non negative real parameter. 
The real function M (oo, t) := lim M(n, t) is shown to admit one unique minimum M(oo, t*) 

rt — >oo 

on the real segment (0,2). Furthermore, the variations of M(oo, t) on the whole real line 
are also studied in detail. 



1 Introduction 

In [3, 4], Itai and Rodeh introduce several symmetry breaking protocols on rings of size N, 
among which the first is considered here. They also show that the average communication cost 
of this particular leader election algorithm takes no more than LN bits, where the value of L is 
computed in [4] to be about 2.441716. 

However, their method is less direct and less general than the asymptotic analysis completed in 
the present paper. Besides, the method is tailor-made for finding only the average number of 
rounds required by the algorithm: the second moment (and a fortiori all other moments), and 
the probability distribution are not considered in [4]. 

By contrast, the asymptotic method used in the analysis of our recurrence relations is very 
general and quite powerful. All moments as well as the probability distribution of the random 
variable can be also mechanically derived from their asymptotic recurrences. A full asymptotic 
expansion, (for large n) can be obtained, and it is illustrated for the mean. An asymptotic 
approximation of the probability distribution (when n — > oo, and j gets large enough) is also 
completed. The latter is derived by computing singular expansions of generating functions around 
their smallest singularity. The present method may serve as a basic brick for finding the com- 
plexity measures of quite a lot of distributed algorithms. 
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The last Section of the paper is generalizing the problem to a probability of the form t/n, 
where t is a non negative real parameter. We show that there exists one unique optimal value 
t* = 1.065439... on the segment (0,2), where the real function M(oo,t) admits one unique 
minimum, M(oo, t*) — 2.434810964 . . ., on the real line. Finally, the variations of M(oo, t) when 
t > 2 are investigated in detail. 



1.1 Algorithm scheme and notation 

For the reader's convenience, we rephrase in our own words the "symmetry breaking" (leader 
election) algorithm designed in [3, 4]. 

Consider a ring (cycle) of N indistinguishable processors, i.e. with no identifiers (the ring is 
said to be "symmetric"), and assume every processor knows N. The leader election algorithm 
works as follows. 

Let n denote the number of active processors. In the first round (initialization), n — N and 
each processor is active. At the beginning of each current round, there remains 1 < n < N active 
processors along the ring. To compute the number of candidates in the round (i.e. all active 
processors that choose to participate in the election), each candidate sends a pebble. This pebble 
is passed around the ring, and every active processor can deduce n by counting the number of 
pebbles which passed through. So, in the beginning of a round every active processor knows n 
and decides with probability 1/n to become a candidate. 

Thus, three cases may happen in a current round: 

• if there is one candidate left, it is the leader; 

• otherwise, the non candidates are rejected (becoming non active), and the remaining active 
processors (the candidates of the current round) proceed to the next round of the algorithm; 

• if no active processors chooses to be a candidate, all active processors start the next round. 

Throughout the paper, we let X(n) denote the random variable (r.v.) that counts the number 
of rounds required to reduce the number of active processors from n to 1 (choose the leader), 
when starting with n~N active processors. The following notations are used. 

P(n,j) := P(x(n)=j), M(n) := e(x(ti)), 

M {2 \n) := E(X(n) 2 ) and <p(n) := E ( e - aX ^ . 

For the sake of simplicity, we also let M(oo) and M^ 2 \oo) denote lim M{n) and lim M^(n) 
(resp.); similarly, P(oo,j) denotes lim P(n,j). 

n^oo 

Finally, let b(n, k) denote the probability that k out of n active processors choose to become 
candidates, each with probability 1/n. In other words, 

\ / 1 



The recurrence equation for the expectation M{n) is easily derived from the algorithm scheme. 
M(n)= 1 + (l-±) n M{n) + £ (") f I )" ( 1 " lT ^ M(k) for n > 1, (1) 



k=2 



and M(l) = (by definition). 
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2 Asymptotic analysis of the recurrence 

Theorem 2.1 The asymptotic average number of rounds required by the algorithm to elect a 
leader is the constant M(oo). When n — > oo, an asymptotic approximation of M(ri) writes 

M(n) ~ - — ^ + Yl ~fcT M W j = 2.441715879 ... (2) 

The second moment of the discrete r.v. X(n) is asymptotically 

M {2 \n)~ - J- 1-1 + 2M(oo) + 5IlT" M(2) ( fc ) ] = 8.794530817..., 



fc>2 



and an asymptotic approximation of its variance (n — * oo) yields 

vax(x(n)j ~ ^ — ^Tj2"( e_1 + (1 - e _1 )5 2 - Slfj = 2.832554383..., 

w/iere Si = V — M(jfe) and 5 2 = V -rr M (2 \k). 

fc>2 k>2 

More generally, 

^»)~ i-e-(° + i) K 1 + E^( fc ) ■ 

y fe>2 y 

Finally, the probability distribution P(oo,j) (n ^ oo) satisfies the following asymptotic approxi- 
mation when j — > oo, 

P(oo,i) ~ r ^z T 2-^ 

w/iere p= .2950911517... 

Up until now, we have been unable to use the classical generating function approach to 
compute M{n). 

However, checking that M(n) is bounded is possible. Indeed, assuming that there exists a 
positive constant B(n— 1) such that 

M(i)<B(n-l) for i = l,...,n- 1, and 5(1) = 0, (3) 

the following inequality holds 

M(n) < —\ -— — ( 1 + B(n - 1) V 6(n, k) ) . 

v ; ~ 1 - (1 - l/n) n - (1/n)™ I ^ i 

So M(n) < B(n), with 



fe=2 
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and -B(l) = 0. (We show below that B(n) is increasing.) 

Let us first analyze the recurrence (4). If B(n) is converging, it must converge to the fixed 
point of Eq. (4), i.e. e. So, we let B(n) = e — A(n), and A(l) = e. 

For fixed k and large n, 




^ - (1 -ir~w. + ^ + ^ + -). «. 



n 



We have 

A(n) = a(n)A(n-l) + (7) 

n 

with 



a(n) 
b(n) 



l-T n -(l/n)"' 
er„_x - 1 

1 - T„ - (l/n)» ' 



Note that n > 3, a(2) = 0, < a(n) < 1/2, and < b(n) < 1. Several constants will be used in 
the sequel: 

e - 2 1 e 1 e - 2 1 e(7e - 13) 

C ° : ~ ' Cl : ~ 2 7^1 ' C2 ~ 2 (e - l) 2 ' ° 3 24 (e - l) 2 ' 

1 -7e 2 + 25e-24 1 c 

C4 := 777 7 Tv? ' c 5 := c i c 2C6 + c 3 , c 6 := , c 7 := — — , c 8 := cic 7 + c 5 c 6 . 

24 (e-1) 2 l-co (l-c ) 2 

For instance, a(n) ~ Co + 0(l/n) and ~ ci + 0(l/n). 
Iterating Eq. (7) gives 

n — 2 n— 2 ,/ .\ z— 1 

A(n) = JJa(n-t)A(t) + E ^ II "(« ~ J") 

i=0 i=0 j=0 

= 1 E TT«(»-i) + E TT«(»-i)- 



Now, 



and so, 



=0 ' j=0 i=n/2 j=0 



E ^ II a (»-i) * 5 E (V2)' - (n^oc), 

»=n/2 j=0 i=n/2 



A(n) - c 6 ci/n. 



Hence, for n sufficiently large, A(n) is decreasing, B(n) is increasing and Eq. (3) holds for n. 
Moreover, A(n) is indeed decreasing to and B(n) converges to e. 
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For the sake of completeness, we can also get a complete characterization of A(n). 
A(n) ~ c A(n - 1) + Cl + = 1} + ^» + ^(« - 1) + /n , 

proceeding by bootstrapping, we first obtain 

oo 

A(»)~ Cl £^ ~ + 

'-^ 7i — % n \ n J 

i=0 

and next, by plugging the above equivalence into Eq. (8), 

A(n) ~ ^ + % + G(l/n 3 ). 

2.1 Asymptotic approximation of M(n) 

Since M{n) is bounded and positive, the limit can be taken in (1) for fixed k, more generally 
for k — o(n 1 / 2 ) (see Subsection 2.2 below). In virtue of Stirling formula and Eqs. (5)-(6), the 
summand writes 

,. e- 1 / k 2 -3k + l 3fc 4 - 22fc 3 + 39k 2 - 9k - 5 \ 
« n '*>~ J^{ l -—2^ + 24^ 

Hence, by Eq. (9), the asymptotic approximation of M(n) is 



M 



\ fe>2 ' / 

which is already given in [4]. 

The average number of rounds required by the algorithm follows, 

M(oo) = lim M(n) = 2.441715878809285246587072 . . . 

n— >oc 

Numerically, 15 terms are enough to obtain a very good precision: the error resulting from the 
sum in Eq. (10) limited to v terms is bounded by 



1-e" 1 ^ k\ 



k>v 

Note also that if the size of the ring is known to be N, the expected bit complexity of the 
algorithm is 2.4417158788 . . . N. It is easily found, since N bits per round are used on the average 
in the algorithm. 

Remark 2.2 Carrying on with the analysis of M(ri) gives mechanically a complete asymptotic 
expansion of M(n). Eqs. (1) and (9) lead to M(n) ~ M(oo) + Ci/n + C^/n 2 + where 

_-i e" 1 ( - k 2 + er Y k 2 + 3fc - 2,e~ x k - 1 + e" 1 - e" 1 ) 

ft - ~w^? + £— w^w> ~ M(k) 



A ,[ + *J + T V - J -M(k) = -.7438715372.. 

4(1 -e- 1 ) 2 ^ 2(1 -e^) 2 k\ 



The expression 0/C2 feeing too long to transcribe, we just give the result: C2 = —.1974635346 . . .. 
The convergence of M(n) to M(oo) is thus very slow: O (n -1 ) • 
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2.2 Interchanging limit and summation 



There remains to justify the interchange of the limit and the summation within the sum in 
Eq. (1), which yields the result in (10). 

2.2.1 Laplace method 

Since the cutoff point in b(n,k) is approximately fco = n 1 / 2 , the asymptotic form of the sum 
^ b(n,k) can be derived from the Laplace method for sums (see [1], [5, p. 130-131]), or 

2<k<n 

"splitting of the sum" technique. 

By taking a suitable positive integer r = o(n 1 / 2 ), we prove that 

n 

i) the sum ^ b(n, k) (the "right tail" of the distribution) is small for large n, and 

k—r 

r i 

ii) lim } b(n,k) — = 0. 

n— *oo ^— ' k\ 

k=2 



i) The ordinary generating function (OGF) of b(n, k), F n (z) := ^ b(n, k)z k is 



fc>0 

1 - z" " 



F n (z) = 

and the OGF of ^ b(n, k) is the product of F n (z) — 1 and l/(z— 1), given by 

r+l<k<n 

F n (z) - 1 
z-l ' 

Considering b(n, k), Cauchy integral formula yields 



r<k<n 



r<k<n 



(z - l)z r 



where is inside the analyticity domain of the integrand and encircles the origin. We see that 
z = 1 is not a singularity for the integrand, so we can neglect the term 1 in the numerator, and 
asymptotically, 



J_ f F n (z) - 1 ^ J_ r exp (n In (l - ^) -rln(z)) 
2th / n Z ~ 2m J n z-l 



dz. 



Again, asymptotically, if we can limit the integration within a neighbourhood of z — 1 = o(n) 
(which is checked below), one obtains 



^ r exp ^ — (1 — z) — rln(z) 
z~l 



27ri X 



dz. 
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To equilibrate, we set z = ry, which yields 



^-t J ^—j exp (ry - r(ha.(y) + In(r))) rdy. 



We now use the Saddle point method. The Saddle point is given by y* = 1 (and z* = r). So we 
set y — 1 + ix and, by standard algebra, we obtain an asymptotic approximation when n — > oo, 



£>(n,fc) 



y2^ r r+l/2 (1 _l /r) ' 

which shows that the right tail of distribution ^ &(n, fc) converges indeed to zero when n — > oo. 
Next, from approximation (9), 



E 



6(n, fc) 



2<k<r 

which tends to zero as n —* oo. 



fc! 



, 2<fe<r 



fc! 



Finally, by completing the sum in (10), it is bounded from above by 



k>r 



k\ 



which also tends to zero as n — > oo. 

Therefore, interchanging the limit and the summation in Eq. (1) is proved justified. 
2.2.2 Lebesgue's dominated convergence method 

The latter justification may also use the Lebesgue's dominated convergence Theorem (see e.g. 
[8, p. 27]). 

By Stirling formula and Eqs. (5)- (6), 



b(n, k) — 



k l,fc 



1 / CX P U ~ 2^ + 2^ 



o 



' n — k s 



Yin) 



k\ k\ 



k\ 



exp 



(i-!r fe+1/2 (i+i^) 



k(k - 3) ^ k l 2k - i - 1 I 



- 1 



k/n 



V 



2n 



2i(i+l) 2n 2n 2 12(1 -k/n) 



(11) 



Set a; = k/n, then 
6(n, fc) 



fc! ~ fc! 



exp 



1 > 
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with 



hix) = (l-x)ln(l-x) + x = y + 0(x 3 ), 

1 3x 
f 2 (x) = - 111(1-*) - x = -— + 0(x 2 ), 

, , . 1- X 1 X 

/s(aO = — s- 



2 12 1-x' 

and 

/i(x)>0, / 2 (x)<0, for|x|<l. 
Thus, for large n, the largest root of nf\{x) + fzix) + fo{x)/n in [0, 1] is given by 

7/n + O (n- 2 ) , 

with 

7 = (3 + V5)/2 = 2.618033988 

which shows that nfi(x) + f2{x) + f^{x)/n > for k > 3 and sufficiently large n (uniformly in 
k). Checking that it remains true for k = n — 5(n), with 5(n) = O (n A ), A < 1, is easy. 

Hence approximation (11) is < for large n, and by Lebesgue's dominated convergence 
Theorem, we can justify the interchange of the limit and the summation in Eq. (1). 

Note that Eqs. (6) and (9) already show that we must take k > 3: the coefficient of 1/n must 
be positive. 

2.3 Asymptotic approximation of M^(n) 

We turn now to the computation of M^(n). 

M( 2 )(l) = 0, and 

M^(n) = Jl--J E((l|I(n)) 2 )| -1 + ^6(n,fc)E((l + A>)) 2 ) 

= 1 + 2 ( 1- ^) M(n) + M(2) W 

n n 

+ 2^6(n,fc)M(fc) + Y, b (n,k)M^(k). 

k=2 k=2 

Hence, when n — > oo (again, interchanging the operators may be justified as in Subsection 2.2), 
M^(n) ~ — L_ f i + 26-^(00) + 2^^M(fc) + ^^M< 2 )(fc) 

e y k>2 ■ fc>2 ■ , 

~ - 1 (-1 + 2Af(oo) + M^ 2 \k)\ = 8.794530817... (12) 



k>2 
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Of course, a full expansion for large n can also be derived step by step. 

Now, since the variance of the r.v. X{n) is denned as var^A(n)^) = M^\n) — ^M(n)j , an 
asymptotic approximation is straightforward (from Eqs. (10) and (12)). 

vaxlX(n)) ~ ~7~ — —r-^ie- 1 + (I - e _1 )S 2 - S?) = 2.832554383..., 



(*(»)) ~ {l _\- 1)2 (e- 1 + (l-^ 1 )^ - S? 

where Si = V — M(fc) and S 2 = V — M {2 \k). 
z — ' k\ z — ' fc 

k>2 k>2 

2.4 Generalization 

More generally, using ip(n) = E (e~ aX< - n ^ as denned in the Introduction, 

<p(n) = e- a ^(l-^j <p(n)+(l-^) ■ 1 + £ b(n, k)<p(k)^ , 



with 



p(l) = 1 and ip(k) = 1 - aM(k) + —M {2 \k) + 



Therefore, 

i-e-(a + i) I 6 '' + 



fe>2 



Also, from the above relations, all moments asymptotic equations can mechanically be found. 

Note that, in contrast to the asymptotic analysis of usual leader election algorithms (e.g. in [2, 
6, 7]), no periodic components are arising in the present asymptotic results. 

3 Asymptotic approximation of P(n,j) 

3.1 Asymptotic recurrence of P(n,j) (n oo) 

The following recurrence on P(n,j) stems from Eq. (1). 

p(„,i) . (i.r 

P(n,j) = 1 - - P(n,j - 1) + ]T &(«, fc)P(fc, j - 1) for j > 1. (13) 



^ \ 71—1 



fe=2 



And the expression of an asymptotic approximation for large n follows, 
P(n,l) ~ e- 1 , 



P(n,j) ~ e-^^j-l) + ^ — P(fc,j-1) forj>l. (14) 

fc>2 
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The above asymptotic approximation on P(n,j) provides the following first 13 values of 
P(oo,j) (j = l, 13): 



.3678794411, .2625161028, .1634224110, .0946536614, .0524658088, .0282518527, .0149122813, 
.0077602315, .0039970064, .0020432067, .0010386252, .0005257697, .0002653262 



Remark 3.1 By definition, the following alternative expressions of M(oo) and M^ 2 \oo) also 
hold, 

M(oo)= ^jP(oo,j) and M^(oo) = ^j 2 P(oo,j). 

So, M(oo) and M' 2 '(oo) comW a/so foe computed from the above definitions. However, more 
than 15 terms should of course be required; viz. about 50 terms are actually needed to obtain the 
same precision as in the previous computations. 

3.2 Asymptotic approximation of P(oo,j) (j — > oo) 

Let us now compute an asymptotic approximation for P(oo, j) when j gets large. First, let 

-l 
k~ 



fe>2 

Whence the recurrence relation (14) also writes 

P(oo,j) = e-^oo.i-l) + D(j-l). 

Here and in the remainder of the paper, the following ordinary generating functions (OGF) 
H(z), G(z) and U(k,z) (of P(oo,j), D{j) and P(k,j), resp.) are used; we define 

H(z) := 2P(oo,jV, G(z) := £l>0V and 

II(fc, z) := P(k,j)z : ', for any fixed integer fc > 2. (15) 

i>i 

From the OGF H(z) defined in (15) and the recurrence (14), we obtain 

H(z) - e~ x z = e- x zE(z) + zG(z), 

and 

z(G(z)+e-i) 
H{Z) ~ l- e -iz ' 
So, i?(-z) has a simple pole at z = e. 

Yet, a numerical check in Eq. (14) shows that P(oo,j) = il (e - - 7 ), and thus, H{z) must have a 
smaller singularity which is (strictly) less than e. 

Now, the OGF II(fc, z) defined in (15) and the recurrence relation (13) yield 
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which gives, for k = 2, 

ikm- 3/2 



l-z/2' 

The above result is of course due to the geometric distribution of P(2,j), with parameter 1/2. 

Hence, 11(2, z) has a singularity at z = 2, and the singular expansion of n(2, z) in a domain 
V around z = 2 stands as 

11(2, z) x l —- . 

V ' 7 l-z/2 

Let R(2) = lim(l — z/2)II(2,z) = 1. In virtue of Eq. (16), it is easily seen that z = 2 is also a 

z— >2 

singularity of all the n(/c, z)'s for any integer k >2. If we denote 

i?(fc) := lim(l - z/2)IL(k,z), 

z—>2 

we derive from Eq. (16) that 

( 1 \ k k 
R(k)= 1-t 2i?(fc) + E 6 (M)2i2(£). 

When k gets sufficiently large, R(k) can be computed (15 terms are quite enough for the precision 
required) . 

Since 

k>2 

the definition of H(z) in (15) shows that z = 2 is also a singularity of G(z). By setting 

e -i e -i 
lim(l - z/2) E n(fc,z) = E'fcT^) =p= .2950911517... 

k>2 ' fc>2 

(again, interchanging the sum and the limit may be justified as in Subsection 2.2), the singular 
expansion of G(z) at z = 2 writes 

Finally, we obtain the singular expansion of H(z) at z = 2, 

{Z >~ (l-2e-i)(l-z/2) ' 

and therefore, when j — > oo, 

P(oo,j) - 2.233499118... 2^. (17) 



4 Numerical results 

As can be seen in the following Figures Fig. 1 and Fig. 2, the previous computations of P(oo,j), 
and M(oo) and M^- 2 \oo) perfectly fit the above ones. Moreover, Fig. 3 shows that the observed 
values of P(oo, j) also perfectly fit the asymptotic approximation of P(oo, j) obtained in (17) for 
sufficiently large j. 
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o : Observed values of P(oo,j) 

— : Asymptotic approximation of P(oo,j) in (17) 



Figure 3: P(oo,j) and its asymptotic approximation in (17) for large j (j — 20, . . . , 30) 



Let t be a non negative real number. Following a question raised by J. Cardinal, let t/n be the 

probability of choosing to participate in the election. 

Is there one unique optimal real positive value t* in some real domain? 

Taken in the initial context of the first leader election ("symmetry breaking") protocol designed 
in [3, 4] (see Subsection 1.1), t is introduced as a real non negative parameter which is assumed 
known to every processor on the ring. 

Initially, all the processors are active. At the beginning of each current round of the election 
algorithm, every active processor knows n (the counting process of n is described in Subsec- 
tion 1.1), and can decide with probability t/n whether to become a candidate in the round. So, 
by definition, t must a priori meet the condition < t/n < 1. 

The recurrence equation for the expectation M(n,t) (with < t < 2) is similar to Eq. (1), 



5 Is 1/n the optimal probability? 




for n > 2 



and M(l,t) = (by definition). 



(18) 
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Upon Differentiating Eq. (18) with respect to t, we obtain 



M'(n,t) = -^1-^ M(n,t)+(l-^ M'(n,t) 



fc=2 



n n \ n 



Now, as in Eq. (10), an asymptotic approximation of M(n, t) for large n yields 

j- A; 

M(oo,i) = 1 + e-*M(oo,i) + ^ e"* — M(fc, i), (20) 

fc> 2 fc ' 



and, similarly, upon differentiating Eq. (20) with respect to t, 

4-fc-l 

M'(oo,t) = -e _t M(oo,t) + e _t M'(oo,t) + ^ e"* — — — M(k,t) 

k>2 

fc>2 fc>2 

or 

M'(oo,t) = 1 - M(oo,t) + e-*M'(oo,i) + ^e~* * M(fc,t) + E e ~* IT M '(M)- 

(21) 

Note that the same expression of M'(oo,t) can also be derived from the recurrence Eq. (19) by 
using asymptotic expansions similar to the ones given in Section 2. 

5.1 Optimal probability on the domain (0, 2) 

A numerical study of the equation M'(oo, t) = on the open segment U = (0, 2) easily leads to 
the solution. 

t* = 1.0654388051 . . . , with M(oo, t*) = 2.4348109638268515517966 . . . 
The relative gain on M(oo, 1) is a bit larger than .0028278945 (hardly more than .28 %). 

Since the (necessary) condition M'(oo, t*) — is not sufficient for M(oo, t) to have an extremum 
at t* , there remains to prove 

1. that M(oo,t) has a minimum at t* G C/, 

2. that this minimal solution t* is indeed unique on the segment (0,2). 
Both results derive from the following Subsection 5.1.1. 
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5.1.1 Af(oo,i) is a strictly convex function on the segment (0,2) 

All definitions regarding real and convex functions that are used in the following may be found 
in [8, Chap. 1 and 3]. 

Since a strictly convex function on some real segment admits at most one global minimum on 
that segment, both above results (1 and 2) are shown simultaneously by proving that M(oo,i) 
is indeed a strictly convex positive real function in U = (0, 2). 

For the sake of simplicity (and in the line of notations in Subsection 1.1), we let M(oo,t) 
denote lim M(n 7 t), 

n—>oc 

(:)(;)' (-r' 



and finally, we also use the notation 

1 



A(n,t) := 1 _ 



for n > 2. 



t/n) n - (t/n) n 

Besides, the following form of the basic recurrence Eq. (18) is considered: 

n-l 

M(n,t) = X(n,t) + X(n,t) ^b(n,k;t) M(k,t) and M(l,t) = 0. (22) 

fc=2 

Starting from the above recurrence Eq. (22), we show below by induction on n, that at any 
point t £ U and for any integer n > 2 all functions M(n,t) are strictly convex positive real 
functions. 

Therefore, as the pointwise limit of such a sequence ( M(n, t) ) in U, M(oo, t) := lim M(n, t) 

\ J n>2 n^oo 

will be itself a strictly convex positive real function in (0, 2) (see [8, p. 73]). 

Note also that, by induction on n, all functions M(n,t) (n > 2) are in C°°(U,M) (i.e., infinitely 
differentiable in (0, 2)), and this is also true for the limit M(oo, t). In the same line of argument, 
M(n,0 + ) := lim M(n,t) = M(n,2~) := lim M(n,t) = +oo for any integer n > 2, which 

remains true in the limit M(oo, t). 

• Basic step. Whenever n = 2, and n = 3, Eq. (18) yields 

,„/„ x 2 n 18-3t-2t 2 

M(2,t)= — r and M(3,i) 



t(2-t) v ' y 3t(2-t)(3-t) ' 

So when fc = 2 and k = 3, M(fc,f) are two positive functions in C°°(U,R) s.t. M(k,0+) = 
M(k,2~) = +oo. 
Moreover, since 

M"(2,t) = 4 3 ^ ( ~ W + 3 4 >M"(2,1) = A and 

„ 2t e + 9t 5 - 189t 4 + 837i 3 - 1674i 2 + 1620i - 648 

M [ ^ t] - ~ 2 3*3(2-t)3(3-t)3 >3 ' 

M(2,t) and M(3,i) are two strictly convex functions in J7. 



15 



• Induction Hypothesis. Assume now that for all t € U, every function M(k,t) is a 
strictly convex positive real function in C°°(U, M), s.t. M(fc,0 + ) = M(k,2~) — +00 for any 
integer 2 < k < n. 

At any point t <G U, A(n, t) > 1 for any positive integer n and 6(n, k; t) > for any pair (k, n) of 
non negative integers. 

71-1 

In virtue of Eq. (22) and the induction hypothesis, A(n, t) b(n, k; t)M(k, t) is a linear 

k=2 

combination of strictly convex (positive real) functions with non negative coefficients, X(n,t) x 
b(n, k; t), in U. 

Furthermore, A(n,t) in infinitely differentiable in U, lim + A(n, t) = +00 and A(n, 2) is bounded 
(except for n = 2, since A(2, 2~) = +00). 

Next, there remains to prove that ^A(n, t)j is also a sequence of strictly convex positive 
real function in U. 

For any given < t < 2 and for any n > 2, the value A(n, t) enjoys the two following inequalities, 
which derive from the tight inequalities shown in [9, p. 242]: for < t/n < 1, 

It is easily seen that, for any fixed value oft G U, A(n, £) (n > 2) is a strictly increasing sequence, 
and lim A(n, t) = . 

rwoo 1 — e~ l 

On the other hand, A(n, t) is a strictly decreasing function of t € U for any fixed n > 2. 

In short, since A"(2,i) > 4 and A"(3,i) > 32/27, A(2,i) and A(3,t) are two strictly convex 
positive real function in C°°(?7, M). 

Again, the proof is by induction on n. If we assume (Induction hypothesis) that, up to any 
integer n > 2, A(n, i) is a strictly convex function of t in J7, then A(n + l,t) is indeed a strictly 
convex function of t on J7. For example, assuming that X"(n,t) > for any integer n > 2, it 
is shown after some algebra that A"(n + l,t) > X"(n,t) > 0, by the above two inequalities in 
Eq. (23) and their resulting properties on A(n, t). 

Thus, the positive sequence I A(n, t) ) is also composed of strictly convex real function in 

V / n>2 

u 

Finally, Eq. (22) and the above results show that, for all t <G U and for any integer n > 2, 
M(n,t) is a linear combination of strictly convex (positive real) functions with non negative 
coefficients: X(n,t) and X(n, t) x b(n,k;t). 

Hence, ^M(n, t)^ is a sequence of strictly convex positive real functions in C°°(U,M), s.t. 
M(n, 0+) = M(n, 2") = +00. 



In conclusion, M(oo, t) is the pointwise limit of the strictly increasing sequence yM(n, t) 

of strictly convex positive real functions of t € U (see [8, p. 73]). Therefore, M(oo, t) is also strictly 
convex in (0,2), and the value M(oo,i*) at t* = 1.065439..., is the unique global minimum of 
M(oo,i) on this segment and we are done. A plot of M(oo,t) is given in Fig. 4. 
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all 



Figure 4: M(oo,i), t € (0,2) 

In that sense, we answered the question set in Section 5: on the real domain (0, 2), t*/n is indeed 
the unique optimal probability for an active processor to choose and participate in the election. 

Remark 5.1 For any integer n > 2, M(n,t) is twice differentiable for all t e (0,2). Hence, if 
M"(n,t) > the functions M(n,t) are all strictly convex; but the converse is not true. 

The positive real function M(oo,t) is defined on the real segment U — (0,2) as the pointwise 
limit of strictly convex positive real functions defined in U . Such is a sufficient condition for 
M(oo,t) to be also strictly convex in U. However the condition is not necessary. 
Furthermore, M(oo,t) is the uniform limit of real functions on any compact subset of the segment 
(0,2). This is another way of deriving that sequences of strictly convex functions do remain 
strictly convex in the limit on any compact subinterval o/(0,2). 

6 What happens to M(oc,t) when t > 2? 

There remains to investigate how M(oo,t) varies as a function of t > 2. In the first place, we 
just assume that the real parameter t belongs to the domain (2, 3). 

6.1 Variation of M(oo,t) in the domain [2,3) 

Since t £ (2,3) and < t/n < 1 (by definition), the value of the function M(n,t) must be 
handled separately in the case when n — 2 (i.e. on a ring with two processors). 
More precisely, two situations may then occur, in which the symmetry cannot be broken with the 
original algorithm (see [4, p. 1]: 

• if t = 2, M(2, 2) = 6(2, 2; 2) = 1. Both active processors on the ring decide with probability 
1 to become candidates in each round, and the protocol either perfoms an election with 
two leaders, or enters an infinite computation; 
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• if 2 < t < 3, we must also set M(2,t) = 1 for the consistency of definitions (when t — > 3~, 
M(2, 3~) = +oo, as is shown below). In such a case no termination of the protocol can be 
achieved. 

Since M(2,t) = 1 is set for all t £ (2,3), the recurrence equation for the expectation M(n,t) is 
expressed in a slightly different form from Eqs. (18) and (22) on the segment [2,3). 

n-1 

M(n,t) = X(n,t) + X(n,t)b(n,2;t) + X(n,t) ^b(n,k;t) M(k,t) and M(2,t) = l, (24) 

k=3 

where, according to the notation in Subsection 5.1.1, 

n\ i t \ I t\ . . 1 



6(n,2;i):= l - 1 , and X(n,t) := -. . , . for n > 3. 

' nj ' v ' ; l-(l-t/n) n -(t/n)" 



There remains to prove that M'(oo,t) > on the segment [2,3), with M(oo,3 ) = +oo. 

First, following Subsection 5.1.1 (i.e. again by induction on n > 3), M(n,t) in Eq. (24) is 
easily shown to be an increasing sequence of n > 3 for fixed t in [2, 3). 
Thus, for all n > 3 and for any t E [2, 3), M(n, t) < M(n, t) < M (oo, t). 

Next, by (modified) Eq. (20) with n > 3 and t € [2,3), upper and lower bounds on M(oo,t) 
are derived. 

More precisely, after few computations the following two inequalities hold for all t G [2, 3), 

M(oo,t) < 26 r + — ^— (25) 
v ; ~ t(t + 2) t + 2 y ' 



M(oo,t) > — -11 + -t 2 e- 1 + M(3,t)e- t {e t -t 2 /2-t- 1) ) 

1 — e ' \ 2 / 



(26) 



9 — 3t 2 — t 3 
where M M = ^ _ f) . 

(Note that since 2.2797 . . . < M (oo, 2) < 2.34726 . . ., both inequalities (25) and (26) make sense.) 

Finally, Eqs. (25) and (26) are used to bound M'(oo,t) from below, and derive that 
M'(oo, t) > on the segment [2, 3). 

Indeed, by (modified) Eq. (21) with t€ [2, 3), a few calculations yield a lower bound on M'(oo, t) 
for any t G [2, 3). 

M'(oo t) > - 2 ^+ t2 ) + 2( e '-t-l)(9-3f-f3) 

M [0 °' t) - tit + 2) t2(t + 2)2 + 3t(t + 2)(3-t) ' 1 j 

And, since M'(oo,t) > 2.26605840... for all f e [2,3), M'(oo,i) > on that segment. 
Furthermore, since all functions M(n,t) (n > 3) are in C°°([2,3),K) (see Subsection 5.1.1), 
M(oo, 3") = +oo holds for all t £ [2, 3). 

Hence, M(oo,t) is strictly increasing on the segment (2,3) and M(oo,3~) = +oo. 
The curve M(oo,t) is represented in Fig. 5 on the segment [2,3). 



18 



Figure 5: M(oo,i), t € [2,3) 



6.2 Variation of M(oo,t) in the domains + 1), with £ > 3 

Investigating the variation of the functions M(oo,i) when i > 3 can be carried out along the 
same lines as in the previous Subsection 6.1. 

As can be noticed (e.g. in Subsection 5.1.1), the only poles of the functions M(n,t) are all 
the non negative integers 0, 2, 3,. . . (1 excepted) on the real line. Thus, the variation of M(oo, t) 
when t > 3 must be considered on all such consecutive real segments (£, £ + 1), where the £'s are 
all integers > 3. 

Since t £ + 1) still meets the condition < t/n < 1 (by definition), each value M(£,i) 
must again be handled separately on each open segment I = (£,£ + !■)■ 

More precisely, whenever n = £ there are £ processors on the ring, and the condition < < 1 
must still hold. The situation is similar to the one in Subsection 6.1: the original algorithm 
cannot break the symmetry, neither if t = £, nor if£<t<£ + l (see [4, p. 1]). 

To overcome the difficulty, and for the sake of the consistency of the definitions, we set 
M(£, t) := |"lg(£)l for all t e /, with £ > 3. For example, M(3, t) := 2 (by definition) on the open 
segment (3,4), and the recurrence for the expectation M(n,t) is slightly different from Eq. (24) 

ifte(3,4). 

Similarly, each basic recurrence equation for M(n,t) (Eq. (24)), M'(n 7 t) (Eq. (19)), M(oo,t) 
(Eq. (20)) and M'(oo,t) (Eq. (21)) must be adapted to the conditions on each segment / con- 
sidered. 

On each open segment J =(£,£ + 1) (£ > 3), the variation of the real function M(oo,i) is 
roughly the same. In particular, M(oo, t) is monotone increasing in /, and it admits no minimum 
on each such segments. 
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7 Conclusions 



As pointed out in the Introduction, performing the asymptotic analyses of various recurrence 
relations brings into play some basic, though powerful, analytic techniques. This is the reason 
why such methods make it possible to find easily all moments of the algorithm asymptotic "cost" 
(the numbers of rounds required), especially M(oo) and M^ 2 \co) (when n gets sufficiently large), 
as well as an asymptotic approximation of P{oo, j) (when j — ► oo). The latter is derived by com- 
puting singular expansions of generating functions around their smallest singularity. Asymptotic 
expansions of all moments can also be mechanically derived. All the numerical results performed 
(with Maple) by both techniques are quite accurate and fit in perfectly. 

Generalizing to a probability t/n, where t is a positive real number, shows that there exists one 
unique minimum of the function M(oo, t) on the real segment (0, 2) : M(oo, t*) = 2.434810964 . . . 
at the point t* — 1.065439 . . . Besides, the variation of M(oo,t) whenever t > 2 shows quite the 
same behaviour on each real open interval (£, £ + 1), where the £'s are all the integers > 2. 

In the asymptotic analysis, the major difficulty arises from the proof of interchanging the 
limits and the summations in the recurrences. Two different methods are given that may be 
used in many other similar situations: the Laplace method for sums, which requires the use of 
asymptotics via the Saddle point technique, and the Lebesgue's dominated convergence property. 

In conclusion, such analytic techniques may serve as basic bricks for finding the asymptotic 
complexity measures of quite a lot of other algorithms, in distributed or sequential settings. 
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